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ABSTRACT 

We present molecular line and 1.4 mm continuum observations towards five 
massive star forming regions at arcsecond resolution using the Submillimeter 
Array (SMA). We find that the warm molecular gas surrounding each HII region 
(as traced by SO2 and OCS) appears to be undergoing bulk rotation. From the 
molecular line emission and thermal component of the continuum emission, we 
independently derived gas masses for each region which are consistent with each 
other. From the free-free component of the continuum emission we estimate the 
minimum stellar mass required to power the HII region and find that this mass, 
when added to the derived gas mass, is a significant fraction of the dynamical 
mass for that region. 

Subject headings: Stars: Formation - Molecular data - ISM: Kinematics and 
Dynamics - sub-millimeter 



1. Introduction 



As a core within an overdense region of a molecular cloud begins to collapse, its angular 
momentum causes the rotating core to speed up, often forming an accretion disk. Until 
recently, it was suggested that for high mass stars, accretion must be halted at the onset 
of an HII region, that the strong out ward pressures must reverse an accretion flow (i.e. 
Garay &: Lizanolll999l : IChurchwelll 120021 ) . However, rec ent models which include the effects 
of gravity have suggested that accretion can continue (jKetol 120031. 120071) and observat ional 
support has followed for infall and accretion in both ionized gas (IKeto fc Woodll2006l ) and 
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molecular gas (jSollins et al.l l2005bl ; iKlaassen fc Wilsonl 120071 . 120081 ) in the presence of HII 
regions. If we are to understand high mass star formation, we must first understand the 
dynamics of their natal environments both before and after the formation of an HII region. 

If analogies can be drawn between low mass and high ma ss star formation, we would 



expec t this continued accretion flow to occur though a disk (i.e. IZhangll2005l ; ICesaroni et al. 



20061 ). If these disks exist, we cannot yet observe them except in the closest sources. However, 
using submillimeter telescopes, we can observe the warm molecular gas around ultracompact 
HII (UCHII) regions to determine whether there is large scale rotation which could translate 
onto smaller scales on which disks may exist. Continuu m studies at high resolution may 



show flattened structures perpendicular to outflows (i.e iGibb fe Hoard 120071 ; IZhang et al. 



20071 ); however, without observing the kinematics of the gas, it is unclear whether these 
structures are indeed rotating. When and how accretion disks form is still fairly uncertain 
for low mass stars, and their presence in high mass star forming regions is still often debated. 
However, there is a growing nu mber of high mass star forming r egions in which disks have 



IRAS 18089 flBeuther fc Walsh 



been detecte d (such as G24A1 (Beltran et al.ll2006T) . IRAS 20126 (ICesaroni et al.ll2005f ). and 



20081 )) all of which have yet to form HII regions. In addition 



to these detections, rotating structu res seen in NH 8 have also been observed in a few sources 
at high resolution with t he VLA fle. lSoUins et alil2005bl : Izhang fc Holll997l : E^eto et al.lll988 



1987 



Zhang et al.lll998l ). Here, we present observations of warm gas tracers (SO2 and OCS) 



towards five high mass star forming regions which are suggestive of bulk rotation in the gas 
surrounding each UCHII region. 

At 220 GHz, the SMA0 in its extended configuration is only sensitive to structures 
smaller than ~ 10". Thus, at the average distance to the massive star forming regions 
studied here (6.4 kpc), we can only observe structures smaller than ~ 64000 AU, and our 
data are ideally suited to probe the warm molecular gas surrounding HII regions. 

In Section [2] below, we br iefly introduce our observations, which come from the same 
dataset as those presented in iKeto et al.l (j2008al ). In Section [3], we present our analysis 
of the velocity gradients in SO2 and OCS and the dynamical masses we derive for each 
region. In Section H] we present our other mass determination methods. In Section 14.11 
we derive gas temperatures based on the J=12-ll K ladder of CH 3 CN. We then use these 
derived temperatures to estimate the gas mass responsible for the SO2 and OCS line emission 
(Section 14. 2p . From the thermal component of the continuum emission we derive a gas mass 
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based on the emission from the warm dust (Section I4.3p . In Section H~4l we present stellar 
mass estimates based on the free-free continuum emission from these sources. In Section 
we discuss the implications of these different mass determinations for each source, and 
conclude. 



Observations 



The observations presented here were obtained at the SMA in its extended configuration 
in the fall of 2005. Source names, pointing centers and source offsets from the pointing 
centers are given in Table [TJ H30ct was centered in chunk 21 of the upper sideband, and 
the observations were made with a spectral resolution of 0.53 km s~ x over the entire 2 GHz 
bandwidth in each sideband. This correlator setup allowed for the simultaneous observation 
of OCS (J=19-18 at 231.061 GHz) in the upper sideband as well as SO2 (J=lli,n-10o,io 
at 221.965 GHz) and CH 3 CN (J=12-ll) in the lower sideband. The K components of the 
CH 3 CN molecule up to K=8 were observed towards each source. CO and 13 CO (J=2-l) were 
also observed; however they are presented in a separate paper (Klaassen et al. in prep). Data 
reduction was completed using MIRIAD in both the uv and image planes. Also presented 
here is the 1.4 mm continuum emission, which was derived from the line free portions of the 
large bandwidth (2GHz) of the lower sideband of these SMA observations. The distances 
to each source and rms noises of these observations are given in Table [TJ The distances 
used here f or the regions in W51 and NGC 7538 IRS 1 are smaller than those used in 



Keto et al.l (12008b), and were taken from the parallax measurements of Ku et al.l (120091 ) and 



Moscadelli et al.l (l2009h respectively. Below, an absolute calibration uncertainty of 20% is 
assumed, but not propagated through our calculations. When uncertainties are listed, they 
arise from the rms in the map, propagated through the equations, unless stated otherwise. 



3. Systematic Velocity Gradients in S0 2 and OCS 

Towards each of our observed HII regions, SO2 and OCS show evidence for a velocity 
gradient across the emission region (see Figures [USD • The solid black lines shown in the 
left hand panels of Figures [IE represent the cuts taken for the position velocity diagrams 
shown in the right hand panels. The directions of the cuts were chosen along the position 
angles (PAs) with the largest velocity gradients, and perpendicular to the outflow direction 
if known. The CO also observed as part of these datasets, when combined with observa- 
tions of the larger scale structures (Klaassen et al. in prep), show outflows from each HII 
region. The directions of these outflows are shown by dashed lines in Figures [U5J The 
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solid ovals presented in Figure [6] show the Full Width at Half Maximum (FWHM) of the 
SO2 emission region deconvolved from the SMA synthesized beam. In both SO2 and OCS, 
the velocity gradients are present in similar size scales and velocity ranges. In two of these 
sources (G10.6 and W51e2, shown in Figures Q] and H] respectively ), si milar velocity gradi- 
ents have previously been observed in NH 3 bv lSollins et all fl2005bh and lZhang fc Hoi (11997V ) 
(respectively). They suggest that rotation is responsible for the velocity gradients, and we 
adopt this interpretation here as well, and extend it to t he rest of our sources. Towards 
G10.6, the velocity gradient seen in NH 3 by ISollins et al.l (j2005bl ) has the same orientation 
and velocity lim its as the observatio ns presented here. The velocity gradient in NH3 towards 
W51e2 seen by IZhang fc Hoi (119971 ) also has the same orientation as that presented here. 



Keto fc Klaassenl (120081 ) see a similarly oriented velocity gradient in the ionized gas (H53a) 



within W51e2 as well. We note that when the dire ction of the bipola r outflow emanating 
fro m these massive sta r form ing regions is known (i.e. iDavis et al.l (119981 ) for NGC 7538 IRS1 
and iKeto fc Klaassenl (120081 ) for W51e2), the orientation of our observed velocity gradient 
is perpendicular to the outflow. The right hand panels of Figures [T]|5] show that the velocity 
gradients in SO2 and OCS. The velocity ranges shown in Figures [T]|5] are much larger than 
the velocity ranges shown in Figure The velocity ranges used for the first moment maps in 
Figure [6] are smaller than those shown in the P V diagrams due to a clipping imposed on the 
data and to highlight the direction of the velocity gradient in each source. The directions 
of the velocity gradients in SO2 and OCS are consistent for a given source, save G28.2, and 
with the exception of W51e8 (see Section |3~2]) . the velocity ranges for both species towards 
a single source are the same. Below, we present our derivations for the velocity gradient in 
each source and our derived dynamical masses. In Section 13.21 we present a caveat to the 
velocity gradient derivation for W51e8. 



3.1. Mass Estimates from Velocity Gradients 



To determine the velocity shifts and angular sizes requi red to calculate our velocity 
gradients, we use the methods presented in lZhang fc Hoi (119971 ). We extracted PV diagrams 
in both SO2 and OCS along the solid lines shown in the left panels of Figures [T]|5j We then 
fit 2D Gaussians to the line emission. In each PV diagram, the Gaussians were fit to the 3 
or 5 a level as stated in the figure captions. The projected major and minor axes of these 
ellipses gave the angular extent and velocity shift across each source. The calculated velocity 
gradients for each source are shown in Table [2] where the uncertainties arise from the velocity 
channel spacing (0.53 km s -1 ) and spatial resolution (0.3") of the observations. 



It is interesting to note that the velocity gradients in G28.2 are differently oriented from 
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each other in the PV diagram, yet result in the same dynamical mass. While the change in 
velocity across the sources is the same for both species, the higher velocity gas is at different 
positions. This is possibly due to the two species being slightly offset from each other as 
seen in the right hand panel of Figure [2J 

If we assume that the material in these regions is in virial equilibrium, we can determine 
the dynamical mass of the enclosed region. Because we cannot constrain the inclination angle 
of these rotating structures from our observations, we will assume that the velocity gradients 
are observed primarily 'edge on' (i.e. the observed velocities are the true velocities of the 
gas) and thus our derived masses are lower limits to the true masses for structures which are 
inclined with respect to the plane of the sky. 

We find velocity gradients in the range of 92-175 km s _1 pc _1 and note that the values 
derived for a given source from S0 2 and OCS agree to within one or two times the quoted 
uncertainties, which is why only the SO2 values are given in Tabled That these values are so 
similar suggest that SO2 and OCS are tracing the same bulk motions. The velocity gradients 
we present in Table [2] are comparable to, although generally slightly smaller than, those 
inferred for three other massive star forming regions (G24A1, IRAS 20126 an d IRAS 18089) 



by app l ying the same analysis meth ods described here t o the data presented by lBeltran et al. 



(120061 ). ICesaroni et al.l (120051 ) and iBeuther fc Walshl (120081 ) (respectively). These sources 



were observed at similar resolutions, but are, in general closer than the sources studied here, 
and do not show evidence for HII regions. We find that our results are generally higher 



than the velocity gradient found in the low mass star forming region L1251 (ILee et al.ll2007l ) 
again using the same analysis on their presented velocity gradients. Our derived velocity 
gradients result in dynamical masses of between 12 and 364 M Q . For all sources except for 
our unresolved source (NGC 7538) the derived dynamical masses are large enough to account 
for one (or more) massive star and surrounding HII region. Because the dynamical mass 
accounts for all of the mass within its outer emitting region (i.e. the hot gas and dust inside 
the HII region as well as the warm molecular gas responsible for the emission), we expect 
these dynamical masses to be the largest values presented here. 



3.2. Outflow Contamination in S0 2 towards W51e8 

In Figure we show both an emission map (left panel) and a PV diagram (right panel) 
for SO2 and OCS in W51e8. In the right panel, the grayscale shows the OCS emission, 
and the contours show the SO2 emission. Here, we see that the velocity gradient in SO2 
around W51e8 extends to higher red shifted velocities than blue (from the rest velocity of 
the source given in Tabled]). The solid ellipse shows the OCS emission. We suggest that 
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the SO2 emission redward of 66 km s _1 is not related to dynamical rotation, and we suggest 
that the SO2 velocity gradient is being contaminated by outflow emission. This excess red 
shifted emission is not observed in OCS (see Panel f of Figure E]); however the orientation 
of the velocity gradient in OCS is consistent with that of the SO2. If we blindly derive 
a dynamical mass for this region from the full velocity gradient in SO2, we estimate an 
enclosed mass of 930 M©. Interestingly the redshifted lobe of the bipolar outflow (as seen 
in CO, Klaassen et al., in prep.) from this source is located near the redshifted part of the 
rotation signature. It is not clear why there is no corresponding CO blueshifted outflow lobe, 
which makes it hard to determine a concrete outflow direction. The correlation in both space 
and velocity between the SO2 emission and CO emission suggests we are seeing outflowing 
material pushing through the SO2 shell of warm gas surrounding the region. The detection 
of H30a in this r egion is marginal a t best, and while it has been detected as a continuum 



ol tldOa m tins r egion is marginal a t best, and wmie it nas been detected as a continuum 
source at 1.3 cm ( Zhang &: Ho! 1997h. it has no t been detected at 3.6 cm (down to an rms 



limit of ~ 1.2 mJy beam" - 1 rtGaume et aDll993h ). These results bring into question whether 
W51e8 is an HII region. Given the large mass of this region however, if it is not an HII 
region, we suggest that it is a very young object which is transitioning from hot core to HII 
region. 

If this high velocity redshifted excess is due to the outflow pushing through the warm 
SO2 emitting region, as suggested by our data, then the molecular outflow observed towards 
this source is continuing to be powered from within the HII region itself. If outflows are 
driven by accretion, this result suggests that accretion is continuing onto the star(s) present 
in this HII region. It is unlikely that thi s red extension in th e emission could come from 



a remnant type outflow (see, for instance iKlaassen et al.l 120061 ) since, by the time the HII 



region would have grown onto observable scales, the outflow mechanism from within would 
have stopped; there would be nothing pushing the S0 2 outwards. We suggest that the 
rotation signatures of the other sources are not contaminated by outflow since the outflows 
from those sources are perpendicular to the velocity gradients present in S0 2 and can be fit 
well with a Keplerian orbit centered on the local standard of rest velocity of the sources with 
no red or blue excesses as seen towards W51e8. 



4. Gas, Dust and Stellar Masses 

In parallel to the dynamical mass estimates from the SO2 and OCS lines, we can estimate 
the mass of molecular gas responsible for the emission (Section 14. 2[) from the integrated 
intensities of the lines. This calculation requires knowing the temperature of the gas in the 
emitting region. The gas temperature for each region was independently derived from the 
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K ladders of the J=12-ll transition of CH3CN. From the K ladders, we also estimate a gas 
mass as described below. 

At 230 GHz, there are two main contributors to the measured continuum emission: a 
thermal dust component and a component from the free-free emission from the HII region. 
In Section 14.31 we discuss how the two components were separated and how a gas mass was 
derived from the thermal dust component. In Section [4741 we discuss how a minimum stellar 
mass for each region was derived from the free-free emission from each source. 



4.1. Temperature estimates from CH 3 CN emission 

CH3CN is a symmetric top molecule, and as such can be used to trace the temperature in 
its emitting region. The individual K levels within a specific J transition (in this case, J=12- 
11) are radiativ ely decoupled, meani ng that the different K levels are populated through 
collisions alone (ISolomon et al.lll97ll ). and the relative populations of the K components 
give an estimate of the ambi ent (kinetic) temperat u re th rough the rotational temperature 
diagram method described in iGoldsmith &: Langerl (119991 ) . We defined the emitting region 
for CH3CN as twice the area of the FWHM of the integrated intensity of the combined 
emission from the K=0-8 components. The CH3CN spectra shown in Figure El are averaged 
over this masked region. In order to determine the integr ated intensity o f each c omponent in 



a syst emat ic way, we followed t he procedures outlined in lRemijan et al.l (120041 ). lAraya et al. 



(120051 ) and lPurcell et al.l (120061 ) which assume local thermodynamic equilibrium and that the 
gas is optically thin. Gaussian profiles were fit to each component, keeping the widths and 
relative rest velocities fixed for a given source. This procedure allows the fitting routines to 
find the best fit to the rest velocity of CH 3 CN in each source. As highlighted in the middle 
and bottom panels of Figure [8], HNCO and CH^CN were within our spectral window but 
not fit by Gaussians. 

From the fitted Gaussians to each K component, we determined the column density of 
the upper level which is related to the integrated intensity of the line through 



N„ 



A u£ hc 3 



Tdv 



where N u is the column density in the upper level of the K component, v is the frequency of 
the transition, A U £ is the Einstein A coefficient for that transition, and J Tdv is the integrated 
intensity in units of K km s _1 . N u is related to the total column density for the emitting 
region through: 
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N u N. 



9u 



TOT c -E u /kT 



Q(T) 



(2) 



where g u is the degeneracy of the upper state, Q(T) is the partition function, A^ TO t is the 
total column density of the spec ies, and E„, is the energy above the ground state of the 
transition (see the appendix of lAraya et al.l 120051 ). Taking the log of both sides of the 
equation gives ln(N u /g u ) oc —E u /kT. The ambient temperature can then be determined 
from the slope of the line of best fit. 

In Figure [7J we show this plot for each of our regions and the derived temperatures are 
shown in Table [31 Since these lines were observed simultaneously, the relative intensities of 
the lines are not affected by calibration uncertainties. The error bars for each data point 
come from the uncertainties in the Gaussian fits, and the uncertainties in the temperatures 
come from the uncertainty in the linear fit to the points in Figure [71 The signal to noise 
ratios of the Gaussian fits to the K=8 components for each source (as well as the K=7 
components for G10.6 and G28.2) were not high enough to be included in the fits, or the 
plots presented in Figure UJ Also, due to the blending of the K=0,1 components they were 
also excluded from the linear fits presented in Figure UJ CH3 13 CN is detected in three of our 
sources. The isotopic ratio of 12 C and 13 C is approximately 55 ([Wilson &: Roodlll994l ). so 
the detection of CHs 1 3 CN suggests our CH 3 CN emission is optically thick. Using Equation 
1 of IChoi et all Jl993h . we find optical depths of <5.8, 50.4, 35.8, <8.9 and 16.3 for G10.6, 
G28.2, NGC 7538, W51e2 and W51e8 respectively. The lower limits on the optical depths 
are for sources in which CH3 13 CN was not detected and the rms noise was used as an upper 
limit on the CH 3 13 CN emission peak. If the gas is optically thick, then we are overestimating 
the temperatures in our regions and underestimating the column density. 

From these data, we find that the temperatures in the warm molecular gas are between 
300 and 460 K. These values are slightly higher than the temperature range in which SO2 is 
expected to be mo st abundant , which is consistent with our derivations giving upper limits to 
the temperatures. iDoty et al.l (120021 ) show that the abundance of SO2 dramatically increases 
at 100 K, and then begins to decline again at temperatures greater than 300 K. This lower 
temperature limit comes from the activation energy required for H2 to combine with OH 
( IDoty et al.l 120021) . At te mperatures greater than ~ 300 K this reaction becomes much less 
efficient ( jCharnleyl Il997l ). The temperature range over wh i ch PC S has a high abundance 



is unclear; however, from the observations of Hatchell et al. 



jl998h_ t appears to have fairly 



similar abundances towards a number of hot cores. iCharnleyl (119971 ) show that after 10 5 
yr, the abundance of OCS should fall as HOCS + is removed from the system (the primary 
source of OCS). 
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The peak velocities of the CH3CN lines are consistent with the Vlsr for each source 
(given in Table [1]), and have similar emitting regions to the SO2 and OCS gas. Thus, 
we suggest the CH 3 CN is tracing the same gas as the S0 2 and OCS, and will use the 
temperatures derived here to estimate the column densities of OCS and S0 2 . 



4.2. Mass Estimates from Molecular line Emission 

We can find, from Equation [2] not only an upper limit to the rotational temperature 
of the CH3CN emitting region, but the column density of CH3CN gas responsible for the 
emission. The slope of the line of best fit gives the t emperature of the region while the 
y-intercept gives ln(N/Q(T)). Using equation A33 from lAraya et al.l (120051 ) we can solve for 
the column density of CH3CN assuming the temperatures derived ab ove for the partition 
functi on. Using the optical depth correction technique described in iGoldsmith fc Langer 
(119991 ) . we calculated the function C r for the optical depths (see their Equations 16 and 22) 
and correct our total column density values accordingly. We use a CH 3 CN abundance of 
2 x 10 -7 (an average of the values presented in lRemiian et al.ll2004l). and we note that other 
authors find lower abundances in Orion (i.e. 10~ 8 — 10~ 9 IWilner et ajj[l994j ), and that these 



abundances would increase our CH 3 CN derived gas masses closer to those p resented below 



from our analysis of SO2 and OCS; however the abundances presented in iRemijan et al. 
were derived based partially on observations of one of the sources observed here (W51e2). 
We present the masses derived from the CH3CN emission both assuming optically thin and 
optically thick emission in Table [3] (the second last, and last columns respectively). 

From the integrated intensity of a single transition of a molecule (such as SO2 or OCS) 
we can also determine the average column densities of the gas responsible for that emission. 
The column density of either SO2 or OCS can be derived by solving Equation [2] for iVroT- 
The column density in the observed species can be converted into a molecular gas column 
density by scaling by the abundance of the species (here S0 2 or OCS) with respect to H 2 . 

For both SO2 and OCS, there are large variations in the abundances given in the 



(Watt k Charnlev 


1985 


Leen & Graf! 


1988; : 


31ake et al. 


1994 


Charnlev 


1997 


Hatchell et al. 


1998|; 


van der Tak et al. 


2003; 


Viti et al. 


200^ 


1) we find S0 2 and OCS abundances of 1 x 10 7 



and 2.5 x 10~ 8 (respectively). The ranges in abundances for these two species were 5xl0~ 8 — 
1 x 10 -6 for SO2 and 1 x 10 -9 — 5 x 10^ 7 for OCS. The mass of each emitting region was 
determined by multiplying the average molecular column density by the area and by the 
mass of molecular hydrogen. We find that, to within our uncertainty estimates, the masses 
derived from the SO2 are approximately twice those derived from the OCS emission. For 
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both species, over estimating the temperature acts to over estimate the mass. We find that 
in general SO2 is observed at a higher signal to noise ratio than OCS, however, given that 



SOo begins to deplete out of the gas phase at temperatures greater than 300 K (IDoty et al. 



20021 ). we may be overestimating the mass from the SO2 emission by overestimating the 



abundance with respect to H 2 . This could be artificially inflating our SO2 derived emission 
masses. These derived masses represent only the warm gas surrounding each HII region, 
and do not reflect the hot material within the HII region, or the cold gas exterior to the hot 
core. Thus, we expect these masses to be lower than the dynamical masses which include 
the stellar and ionized gas masses within the HII region. 



4.3. Mass Estimate from Continuum Emission 



The two main components of the continuum emission in our sources are free- free emission 
from the HII region, and the warm dust from its surroundings. The predicted emission 
from the free-free component of the continuum flux can be calculated via equation 2.124 of 
Gordon Sorochenkol (120021 ) . This free-free emission can then be removed from our observed 
1.4 mm continuum, and we assume that the remaining continuum emission is due to thermal 
dust. 



Using the equations and relations presented in iHildebrandl (119831 ). we determined the 
dust masses in the region s producing the d ust continuum emission. We have extrapolated 
the dust emissivity from IHildebrandl fjl983h as Q x = 7 x 10- 4 (0.125/A OTm )^, where 0.125 
and X m m are both wavelengths in units of mm (here, A mm = 1.4). This expression was 
derived without assuming that (3 = 1 at 250 /jm. For our calculations, we assume a gas to 
dust ratio of 100 which allows for the conversion from dust mass to gas mass. We further 
assume that the gas and dust temperatures are not well coupled within the HII region, since 
the gas temperature within the HII region is well above the dust sublimation temperature. 
Instead, we assume tha t the central star is heating the dust, and the temperature falls off as 
Tdust = {R*/R) 0A x (ILamers fc Cassinellilll999l ). For the region in which the temperature 
is below the dust sublimation threshold, the temperature gradient was then averaged over 
our 1" beam (or, a linear size of 6400 AU at the average distance to our sources) to give an 
average temperature of 400 K. The masses listed in Table [2] were calculated using this single 
temperature and assuming (3 = 1.5. Increasing (3 to 2 (while keeping the temperature at 400 
K) increases the derived masses by more than a factor of 3. If we use the gas temperatures 
derived from our CH3CN the changes in the masses derived from the dust emission are within 
one sigma of the values listed in Table [2j 



We find, in general, that the dust masses derived in this manner (and presented in Table 
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[2]) are lower than the derived dynamical masses presented in Section 13.11 In fact, we find 
that the average dust derived mass is approximately half of the average dynamically derived 
mass. This result is expected since our calculated dust masses are for regions which are not 
quite as extended as the SO2 emitting region and do not reflect the masses of the central 
star(s) as is the case for the dynamical masses. 



4.4. Stellar Mass Estimates 



Using the spectral energy distribution fits to the continuum data presented in lKeto et al.l 



( I2008al ) (both with and without density gradients within the HII region) we estimate the 
masses of the central stars powering the observed HII regions. From their continuum mea- 
surements, the ionizing flux required to maintain the HII region was determined and a 
minimum stellar mass inferred. If the HII region were expanding, the stellar mass required 
to produce the necessary ionizing flux would be higher. The required minimum stellar masses 
to maintain these HII regions (assuming density gradients within the HII region) are pre- 
sented in T able [2] for all sources except W51e8 which was not included in the analysis of 



Keto et al. 



As stated above, we find that the average source mass derived from the continuum 
emission is approximately half of the dynamically derived mass. When the stellar mass is 
added to the dust mass, we find that the combined mass is within approximately 1 a of 
the dynamical mass for G10.6 and G28.2. For NGC 7538, the combination of dust and 
stellar mass is much greater than the dynamical mass; hence the estimate of 12 M Q for the 
dynamical mass must be an underestimate. For W51e2, we find that the dynamically derived 
mass is much greater than the combined dust and stellar mass estimates, which suggests that 
perhaps the dust derived mass maybe underestimated. 



5. Discussion and Conclusions 



Using SO2 and OCS emission, we have detected rotational signatures in the warm 
molecular gas surrounding five HII regions. The bulk rotation of the gas inside the HII 
region (as seen within W51e2 by lKeto fc Klaassenll2008l ) and around them (as shown here) 
suggests that each of t he HII regions may have formed from a single collapsing core (i.e. 
Reid &: Matthews! 120081 ). If these rotational motions are indicative of a smaller scale disk 
within each HII region, these disks m ay be simila r to the photoevaporating disks within 
hypercompact HII regions described in iKetol (120071 ). 
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Here we have presented two methods for determining the gas masses of the hot cores sur- 
rounding each HII region: from the molecular line emission and from the thermal component 
of the continuum emission. Generally, the SO2 emission mass is larger than the dynamical 
mass, however, as stated earlier, we may be over estimating the SO2 abundance at the high 
temperatures we derive for our sources (which are upper limits due to optical depth effects), 
and our dynamical masses are lower limits since we do not know the inclination angles of 
the individual sources. We find that the average of the derived warm gas masses (S0 2 and 
OCS for a given source), when combined with the average minimum stellar mass required 
to power the HII region, are a significant fraction of the average derived dynamical masses. 
Note that for NGC 7538, the stellar and gas masses are much larger than the dynamical 
mass, suggesting the dynamical mass is underestimated for this unresolved source. That 
the gas and stellar masses are a large fraction of the dynamical mass in each resolved region 
suggests that these regions may be in rough dynamical equilibrium. While we find that most 
of the continuum emission in each source is due to free-free emission, it is interesting to note 
that the expected dust contribution to the continuum emission is the highest in W51e8, a 
region we suggest to be quite young and only beginning to form its HII region. In our other 
four sources, the dust emission contributes less than half of the continuum emission, while 
in this source it contributes at least 72% of the emission. Othe r authors have also suggested 
that this region may be a mostly dusty source (jLai et al.ll200ll ). 



NGC 7538 is the only source in this sample in which the deconvolved SO2 emitting 
region and continuum source are not resolved by the SMA beam. It is also the only source 
where the dynamical and gas masses are too small to be responsible for forming the massive 
star(s) sugges ted by the stella r mass es presented in Table El The outer radius of the maser 
disk seen by iPestalozzi et al.l (120041 ) in this region is ~ 1000 AU, also smaller than our 
resolution element. That the emission is unresolved greatly increases the uncertainty in our 
mass estimates. 

We note that previous s tudies have found m uch higher masses for the G10.6 system (i.e. 
a few hundred solar masses, iKeto &; Woodll2006l ). however the masses we derive here appear 
to be internally self consistent. In contrast, the cores in W51e appear to be much more 
massive than G10.6. This may be due to the active hot cores surrounding the HII regions 
in these sources enhancing both the chemical abundances and radii of the sources; each of 
which act to increase our mass estimates. 

We find the masses de rived for the two reg ions in W51e to be quite large. The dynamical 
mass of W51e2 derived by IZhang fc Hoi (119971 ) is 10 M on size scales of 0.3". Our value of 
364 M Q within 1.2" is consistent with their mass projected out to the size scales we observed 
(M2 = Mi * (V^^/V^ri)), given the uncertainties in the measurements. The large masses 
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associated with these two regions suggests that they must be forming massive clusters of 
stars, as has suggested to be the case for G10.6 as well (i.e. ISollins et al.ll2005bl ). 

As discussed earlier, there appears to be a component of the SO2 emission in the W51e8 
region which may be due to outflow. That we detect outflowing gas in the warm layers 
surrounding the HII region suggests that the mechanism driving the outflow must be enclosed 
within the shell of SO2 emission; it likely originates somewhere inside the HII region. This 
supports the suggestion that there is ongoing accretion onto this (these) massive protostar(s). 



We would like to acknowledge the support of the National Science and Engineering 
Research Council of Canada (NSERC). We would also like to thank our anonymous referee 
whos insights greatly improved this manuscript. 
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Table 1. Observational Properties 



Source 




Coordinates 21 


VLSR b 


Dist. 


rms c 


Beam 


Source Size d 




RA 
(hms) 


DEC Aa 

r") n 


AS 

(") (km s- 1 ) 


(kpc) 


(Jy/Beam) 


maj.xmin. PA 

(") CO (°) 


maj.xmin. PA 

(") (") (°) 



G10.6-0.4 18 10 28.7 -19 55 49 11 15 -3 

G28.2-0.04 18 42 58.17 -04 13 57.0 -1.5 -0.3 97 

NGC 7538 19 13 22.069 10 50 52.5 0.3 0.3 -59 

W51e2 19 23 43.913 14 30 14.7 0.6 -0.3 55 

W51e8 19 23 43.913 14 30 14.7 -6.6 55 



6.0 


0.32 


1.5x0.9 


61.0 


4.8x1.8 


-64.9 


9.0 


0.19 


1.3x1.0 


-84.5 


1.5x1.1 


-16.0 


2.7 


0.12 


1.2x0.9 


-79.3 


0.9x<0.9 


15.8 


5.1 


0.23 


1.3x0.7 


-88.6 


2.8x1.9 


17.5 


5.1 


0.23 


1.3x0.7 


-88.6 


3.1x2.6 


19.5 



Coordinates given are for the pointing center of the observations; the offsets show the source center with respect to i^je 
pointing center. 



b V T ,sR values tak en from: IPurcell et all (120061 ) (G10.6), iSollins et all (l2005ah (G28). Izheng et all feOQlh (NGC 7538 IRS1), 



Sollins et all (12004 ) (both W51e regions). 



C RMS noise limits refer to channel spacings of 0.53 km s 1 at 220 GHz. 
d In SO2 emission deconvolved from synthesized beam. 



Table 2. Dynamical Masses, Dust Derived Masses and Minimum Stellar Masses 



S0 2 Dust derived Masses Stellar Masses 



Source 


P.A a 

(°) 


Radius b 
(") 


Vel. Gradient 

(km s _1 pc _1 ) 


Mass 
(M©) 


odust 
° 1.4mm 

(Jy) 


Mass 
(M©) 


(M ) 


G10.6-0.4 


140 


1.5 


92±15 


162±27 


2.0±0.2 


136±14 


69 


G28.2-0.04 


173 


0.45 


153±26 


68±20 c 


0.6±0.2 


95±34 


34 


NGC 7538 


223 


1 


156±44 


12±3 


2.5±0.7 


36±10 


62 


W51e2 


51 


1.5 


175±22 


364±51 


3.9±1.5 


140±53 


87 


W51e8 





1.6 


157±25 


320±45 


1.8±0.7 


82±33 





a Position angle of the cut for the PV diagram. 

b Estimated radius of the rotating region. Uncertainties are estimated to be half a pixel 
(0.15"). 

c Dynamical mass estimate for this source using OCS emission is 41±17 M . 



Note. — isted uncertainties are due to measurement error, not calibration uncertainty, which 
we estimate at 20% 



Table 3. Temperatures and Molecular Emission Derived Gas Masses 



ooui ce 


T, a 
J- kin 

(K) 


\/p1 T-? ^ n crp c 

V CI. J-LfXllciC 

(km s _1 ) 


f Tdv 

(K km s- 1 ) 


(10 15 cm" 2 ) 


Mxi b 
ivih 2 

(M ) 


ivih 2 

(M©) 


lnC NZO(T)) 


PU PAT 

iN CH3CN 

(10 15 cm~ 2 ) 


Mu b 
ivih 2 

(M ) 


M r>1 

iVi H 2 

(M ) 


G10.6-0.4 


323±105 


-6,4 


113±15 


37±4 


85±11 


57±3 


24.2±0.2 


1±0.3 


1±0.3 


5.6±1.6 


G28.2-0.04 


300±89 


90,100 


175±27 


52±7 


131±20 


65±6 


24.4±0.1 


1±0.7 


1±0.6 


47±28 


NGC 7538 


310±94 


-65,-45 


297±12 


92±3 


7.6±0.6 


3.1±0.4 


24.2±0.2 


1±0.6 


0.2±0.1 


5.9±3.6 


W51e2 


460±130 


45,75 


799±58 


429±27 


525±45 


294±32 


26.3±0.1 


21±14 


16±11 


265±188 


W51e8 


350±83 


45,75 


634±58 


232±18 


385±37 


220±27 


25.9±0.1 


8±4 


7±4 


66±34 



I 



a Upper limit on gas temperature due to optical depth effects in the CH 3 CN emission ^ 
b S0 2 , OCS, and CH 3 CN abundances used to calculate the total gas mass were lxlO- 7 , 2.5xl0~ 8 , and 2xl0- 7 (respectively) 
c isted velocity ranges were used to create the zeroth and first moment maps 



Note. - 



isted uncertainties are due to measurement error, not calibration uncertainty, which we estimate at 20% 
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Fig. 1. — (G10.6-0.4) Left: 1.4 mm continuum emission cut at 3a is shown in grayscale with 
SO2 and OCS emission shown with solid and dotted contours (respectively). The solid line 
shows the cut used for the PV diagram shown on the right while the dashed line shows the 
direction of the outflow as observed in CO (Klaassen et al. in prep). Contours begin at 3a 
for each integrated intensity map, and increase in increments of 2a. Right: Position- Velocity 
(PV) diagram of OCS (greyscale) and SO2 (contours) emission along the cut shown in the 
left panel. Again, the greyscale and contours both start at 3a, and the contours increase in 
intervals of 2a. The solid oval shows the fit used to derive the velocity gradient in the source 
from the S0 2 emission. 
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Fig. 2. — (G28.2-0.04) This figure is similar to Figured! however the contours in the left 
panel start at 5a instead of 3. The two ovals in the right panel show the fits to the SO2 
emission and the OCS emission. This is the only source for which the two species do not 
show the same velocity gradient. 
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Fig. 3. — (NGC 7538 IRS1) This figure is similar to Figure [IJ however the contours in 
the left panel start at 10cr instead of 3, and increase in intervals of 10a instead of 2. In the 
right panel, the contours begin at 5a and increase in intervals of 5a. 
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Fig. 4. — (W51e2) This figure is similar to Figured], however the contours in the left panel 
start at 5a instead of 3, and increase in intervals of 4cr instead of 2. 




Fig. 5. — (W51e8) This figure is similar to Figured], however the contours in the left panel 
start at 5er instead of 3, increase in intervals of 4cr instead of 2, and the oval represents the 
OCS emission instead of the S0 2 emission. Note that only the red outflow lobe has been 
seen in CO for this source (Klaassen et al. in prep). 
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Fig. 6. — First moment maps of SCV Lines show the direction of the velocity gradient, 
while thick ovals show the deconvolved SO2 emitting regions (FWHM), and thin contours 
show the continuum emission starting at 5cr and increasing in intervals of 5a (except for 
NGC 7538 where the contours start at, and are in intervals of 10cr). Panels e and / show 
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Fig. 7. — Plot of the natural logarithm of the column density in the J=12-ll transition of 
CH3CN divided by the degeneracy of the state as a function of the energy of the observed K 
component. The inverse of the slope of the line of best fit shown here gives a measure of the 
temperature in each region. The K=0,1 lines were blended in each spectrum, and thus the 
column densities for these components were not included in the fits. For most sources, the 
integrated intensity of the K=8 component was below the noise level, and thus these values 
were also not included in the fits. 
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Fig. 8. — CH 3 CN integrated intensities and gaussian line fits for each source. The vertical 
lines denote the expected velocities of each K component, and the velocities plotted along 
the x axis are with respect to the local standard of rest velocity. The plots have been aligned 
according to the rest velocity of the sources listed in Tabled! The HNCO and CHg 3 CN lines 
near the K=6 component of the CH 3 CN line have been marked in the middle and bottom 
panels to show that they are contributing little contamination to the integrated intensity of 
the K=6 transition. Note that in NGC 7538, an inverse P-cygni profile can be seen in each 
K component suggesting that the CH3CN may be tracing an infalling component of the gas. 
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